Modeling of microwave-assisted switching in micron-sized magnetic ellipsoids 
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The microwave assisted magnetisation reversal is modelled in a permalloy micron-sized magnetic 
ellipsoid. Our simulations confirm that this process requires less field than magnetisation reversal 
under a static field. This is due to a different reversal mode which in case of the microwave-assisted 
process is always a ripple structure. During the magnetisation reversal two stages: nucleation and 
relaxation are distinguished. The nucleation process is governed by spinwave instabilities. The 
relaxation process is related to the domain expansion through domain wall propagation determined 
by the precessional motion of magnetic moments in the center of the domain walls. As a consequence, 
the switching time is a complex oscillating function of the microwave frequency. 

PACS numbers: 75.40.Gb 75.40.Mg 75.75.-l-a 



I. INTRODUCTION 

The fundamentals of fast magnetisation switching con- 
stitute one of the important problems in magnetism, 
related to technological applicationsi. Experimentally 
microwave-assisted switching has been reported as re- 
quiring a smaller applied field thus opening new possi- 
bilities for applications in magnetic recording and com- 
munication devices^i^i^'^i^. The microwave-assisted (mw) 
switching process with linearly-polarized excitation can 
be realized based, for example, on a fast magneto-optic 
Kerr effect set-up, applying simultaneously static and mi- 
crowave frelds^i^i^ii. Another possibility to observe mw- 
assisted switching is given by measuring the magnetore- 
sistance as has been done, for example, in Co stripes^ or 
in NiFe magnetic tunnel junctions^. These experiments 
use relatively large (up to tens of microns) lithographi- 
cally prepared magnetic elements whose dimensions do 
not allow the occurrence of homogeneous magnetisation 
processes. 

Differently from this, the micro-SQUID experiment*^, 
performed on Co nanoparticles with 20 nm diameter 
gives an example of a macrospin magnetisation be- 
havior. In this system nonlinear effects are clearly 
observed experimentally in agreement with a simple 
macrospin dynamical model. The behavior found in 
the macrospin magnetisation dynamics has stimulated 
the technologically important proposal of micro-wave as- 
sisted (or ferromagnetic resonance (FMR) assisted) mag- 
netic recording^iiSiii. 

The process of understanding the underlying physics 
is still going on and several ideas arc found in the liter- 
ature. From the theoretical point of view, probably the 
most comprehensive understanding has been achieved in 
the case of one single magnetic momenl^^ d'-'i^^d^d'^d^d^ . 
In what follows we would like to mention some of the 
proposed mechanisms, since most of them can also be 
partially encountered in the complex situation of non- 



homogeneous magnetisation reversal processes discussed 
in the present article. The most frequently used one is 
based on a simple idea of the efficient energy transfer dur- 
ing excitation with frequencies close to the FMR one^i^^. 
It has been also argued that excitation with a circular 
polarised field would be much more efficient due to the 
possibility to synchronize with the rotation sense of the 
magnetic momen1>i^. Although it seems to be less ef- 
ficient, synchronization is also possible with a linearly 
polarised fieldi^. This idea remains certainly valid in a 
general case, however it should be noted that the notion 
of the FMR frequency, as corresponding to excitation of 
mode with zero wave vector, is not exact for micron- 
sized magnetic elements. Due to the minimization of 
the magnetostatic charges at the surfaces, even small an- 
gle magnetisation vibrations are not homogcneouSi^ and 
the frequency spectrum depends strongly on the size and 
shape of the magnetic nanoelementi^iil. Moreover, in 
magnetic nanoelements the frequency spectrum becomes 
quantize d^^d^ . The second idea found in the literature 
is also simple: The microwave field provides a constant 
energy input which rises the magnetic energy of the par- 
ticle up to a level higher than the saddle point of the 
magnetic energy landscap e^i^ i 1 2 , 1 5 _ Note that here also 
the thermal-activated process becomes possibl o^'^i^° . In 
this sense, it can be said that the microwave field can 
act as an energy source which effectively decreases the 
energy barrier. It should be noticed here that this idea is 
relied on a simple picture of a two-level system which is 
generally not valid for a non-homogeneous magnetisation 
reversal scenario. 

The fast mw-assisted magnetisation switching pro- 
cess in magnetic elements is closely related to the phe- 
nomenon of precessional switching. Indeed, a small per- 
pendicular field helps the switching to take place via the 
excitation of the precessional motion2ii22ii^ due to the 
torque acting on the magnetisation. The main problem 
is that with, for example, a perpendicularly applied field 
the magnetisation continues to precess ("ringing" phe- 
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nomenon) and could switch back^l. Furthermore, it has 
been suggested to use a pulsed field and to tune the pulse 
duration to avoid multiple switching. Even the pulse 
shape could be optimized^, together with the use of two 
pulses with an optimised delay time^ii^. This case again 
is mostly studied in the case of one magnetic moment, al- 
though micromagnetic simulations have been also widely 
performed in nanoelements with different shapes^l. 

The final idea which we would like to mention is that 
one magnetic moment under a microwave field represents 
a classical example of a parametrically excited nonlinear 
oscillator. Thus the problem is similar to the classical 
parametrically excited pendulum, where parametric in- 
stabilities occur varying the strength and the frequency 
of external oscillation. Such instabilities may be con- 
sidered as precursors of the switching process. Next, in 
these systems nonlinear phenomena such as bifurcations 
and chaosi^i^i2^ occur. A classical example of this is 
the appearance of an additional large-amplitude station- 
ary orbit via the folding bifurcation, again studied in the 
case of one macrospin onl y^^'^^ . The nonlinear phenom- 
ena lead to a complex structure of the trajectories near 
the separatrix, a special trajectory separating the basis of 
attractions of the two minima. Under the influence of ex- 
ternal periodic force, the separatrix may become fractal 
so that closely situated initial conditions may lead either 
to switching or not^^. This produced, for example, fractal 
regions of switching beyond the Stoner-Wolfarth astroid 
observed experimentally in micro-SQUID experiments^ 
and numer ically in Refsi[llll[ll. However, since these 
phenomena are complex even in the case of one magnetic 
moment, in a large system nonlinear phenomena are more 
difficult to interpret and may not lead to a clear pattern. 

The reversal under oscillating fields of large-sized 
magnetic nanoelements was mostly studied in cases 
with simple magnetisation reversal modes, such as do- 
main wal&^, C or 5 magnetisation states^, the vor- 
tex to onion state in magnetic rings^ or vortex core 
reversa&^. The enhanced mobility of domain walls 
with some frequencies was numerically found^. Exper- 
imentally, using the magnetic circular dichroism tech- 
nique in a magnetic element having a Landau pattern, the 
asymmetry of the domain wall motion under a linearly- 
polarised mw field with the drift in one direction was 
reported^^. This has been explained by the entropy max- 
imization involving the precessional excitation (the exci- 
tation of a so-called " self-trapping spin- wave mode" ) . 

As we mentioned above, most of the previous theoreti- 
cal studies involve either fast switching in the macrospin 
case (or effectively conditions of coherent rotation) or 
the switching of individual objects as vortices or domain 
wall structures. However, the spatial resolution of the 
fast Kerr experiments^ii^ requires quite large magnetic 
elements with not so simple and unique magnetisation 
reversal modes. Namely, the mw-assisted switching in 
large systems produces a sequence of nucleation - prop- 
agation - relaxation phenomena which is the subject of 
the present study. We show that the overall dynamics is 



a complex phenomenon where some of the above men- 
tioned processes could occur simultaneously and could 
play a role on different stages of the magnetisation rever- 
sal. 



II. MODEL 

In the present work and with the aim to understand 
the mw-assisted switching processes in large magnetic 
elements^ii^, we use a micromagnetic model to simulate 
the magnetisation dynamics in a permalloy ellipsoid. Due 
to the limitation of the simulational size we used reduced 
dimensions of the ellipsoid: the total magnetic volume 
is y = X 2/xm by 25nm. It has the same as- 

pect ratio as in the experiment. We used the publicly 
available code Magpar— with a finite element discreti- 
sation and in which we implemented the possibility to 
simultaneously apply dc and ac fields. The permalloy 
easy axis was directed parallel to the large ellipsoid axis 
(y-axis, see Fig. [TJ. Initially the ellipsoid was magne- 
tised along this axis. The dc-applied field is placed par- 
allel to it in the opposite direction. The linearly po- 
larized ac field with large amplitude is directed perpen- 
dicular to it (x-axis). The following micromagnetic pa- 
rameters were used in the simulations: anisotropy value 
K = 3.31 X lO^erg/cm^ , saturation magnetisation value 
Ms = 860 emu /cm? (the anisotropy field Hk = 7.70e), 
exchange parameter A = 1.05 x 10~^erg/cm, damping 
parameter a = 0.012. 
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FIG. 1: Geometry of the simulated system. 

In Fig. [5]we present a hysteresis cycle for the permalloy 
ellipsoidal element with the field applied parallel to the 
easy axis direction. The coercive field is much larger than 
the anisotropy field value due to the large magnetostatic 
shape anisotropy. 
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FIG. 2: Simulated descending branch of the hysteresis cycle 
for the permalloy ellipsoid with static applied field only (full 
circles) and with simultaneously applied static and mw field 
with Hac ~ 25.10e, and two frequencies (open symbols). The 
dc-field is applied along the long ellipsoid axis Y. The average 



ellipsoid magnetisation 
saturation value. 



< Aly > is normalized to the total 



are almost independent of the ac-field frequency which 
role is just to induce the precession. The results arc dif- 
ferent when strong deviations of the magnetisation occur 
which stresses the importance of the nonlinear effects. 
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III. RESULTS 



FIG. 3: Amplitude of the critical field necessary for mag- 
netisation switching as a function of the maximum field an- 
gle with the anisotropy axis for static field and microwave- 
assisted switching. 



First, wc model the hysteresis cycle with simultane- 
ously applied dc and ac fields (see Fig. ^ and note the 
reduction of coercivity in the case of the mw-assistcd 
switching process. The hysteresis cycle is square-like in 
the former case and rounded in the latter one. To un- 
derstand the coercivity reduction we model the minimum 
strength of the applied field required to produce the mag- 
netisation reversal in our system for different mw-field 
frequencies. 

We notice that one of the large contributions to the co- 
ercivity reduction comes from the deviation of the mag- 
netisation direction from the y-axis, following the ac- 
field. Thus, the coercivity is reduced simply due to the 
fact that most of time the resulting field is applied at 
some angle to the easy axis. Consequently, we com- 
pare the mw-assisted switching process with the situa- 
tion when a static field is applied at an angle equal to 
that formed during the mw-assisted case with the maxi- 
mum amplitude Hac of the ac field. The maximum am- 
plitude of the critical field Her = a/ H^^ + H^^, necessary 
to switch the magnetisation, is presented in Fig. [3] as a 
function of the maximum applied field angle. We use 
the following condition for the magnetisation switching 
< My >< 0, where < My > is the average magnetisation 
along the y direction. Our results show that even with 
this more fair comparison, mw-assisted switching always 
requires less field, except for the case of very large devi- 
ation angles and some frequencies. Therefore, the first 
important contribution to decrease the coercivity dur- 
ing the mw-assisted switching process is the deviation 
of field from the easy axis which leads to precessional 
switching. Moreover, for small field angles the results 



To illustrate the differences between the static and 
mw-assisted switching processes we present the dynam- 
ical configurations during the magnetisation reversal in 
Figs, m and [5l For the static process (Fig. |4|) the mag- 
netisation nucleation starts with two vortices in the op- 
posite upper and lower sides of the ellipsoid. During the 
irreversible jump these vortices propagate in opposite di- 
rections creating two domain walls. The propagation of 
these domain walls toward the left and the right sides of 
the ellipsoid completes the magnetisation reversal. 
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FIG. 4: Magnetisation configurations during the static field 
hysteresis process at different time moments (a) initial state at 
the coercive field, (b) intermediate state, showing the vortices 
propagation and (c)the final state at the coercive field. The 
field is applied at 45° with the Y axis. 

Fig. [5] shows the magnetisation configurations during 
the mw-assisted switching process. The complimentary 
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video material (see Ref. IstI ) illustrates the temporary dy- 
namical evolution of the demagnetisation reversal. We 
observe that the ellipsoid is dynamically divided into 
domain structures. The number of domains depends 
strongly on the mw frequency and may be related to the 
length of the spinwavc mode, excited in the system. In 
Fig. [S]we represent an approximate number of magnetisa- 
tion domains nucleated during the first two nanoseconds 
as a function of the ac-field frequency and in the mag- 
netisation range 0.6 < < My > < 0.8 (normalized to the 
total saturation magnetisation). Note that the domain 
size is larger in the central part of the ellipsoid than in 
the upper and lower sides, due to the ellipsoidal form, 
and initially grows with time (see Fig. [7]). The division 
of the ellipsoid into domains during mw-assistcd switch- 
ing is experimentally confirmed by Kerr images in Ref.l 
Therefore, the mw-assisted case is characterized by a re- 
versal mode completely different from the one appearing 
during the static hysteresis. 




FIG. 5: Dynamical configurations: (upper) and Mz 

(lower) components (grey scale) during the mw-assisted mag- 
netisation reversal for applied fields Hdc = — 31.410e, Hac ~ 
18.840e and for two different frequencies (left) v — IGHz and 
(right) u — QGHz. For the dynamics of the reversal process, 
see the supplemental material^S. 

Generally speaking, the ac-driven process is compli- 
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FIG. 6: Number of domains in the ellipsoid nucleated during 
the first two nanoseconds as a function of mw frequency for 
applied field H^^ = -31.410e, Hac = 25.130e. 
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FIG. 7: Domain size in the ellipsoid along the y coordinate 
as a function of the position of its center for two time snap- 
shots, corresponding to the first nucleation-expansion stage of 
the mw-assisted reversal processes and for Hdc ~ —31.4106, 
Hac ~ 25.130e. Note that after approximately t > 2.5ns the 
domain size expansion is suppressed. 



cated and involves different objects. The domains are 
separated by domain walls (of Neel or cross-tie types) 
and the magnetic moments in the center of these domain 
walls are constantly precessing (see Fig. [5] (c) and (d)). 
This precession is important for the domain wall mobil- 
ity. In the junctions between the domains vortices are 
created (see Fig. [5]). Next, during the process we ob- 
serve constant spinwave generation and their reflection 
from the ellipsoid boundaries and domain walls. There- 
fore, the overall process contains an interplay of several 
important effects, discussed earlier in the literature in 
more simple situations (see the Introduction). Each of 
them has been reported previously to be influenced by 
the mw-field. 

Fig. [5] shows the dependence of the magnetisation 
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FIG. 8: More detailed magnetisation configuration in tlie 
upper part of tfie ellipsoid showing the occurrence of do- 
main walls and vortices. Hdc ~ — 31.410e, Hac = 25.130e, 
u = O.lGHz. 



switching time in the ellipsoid on the mw frequency for 
several amplitudes of the ac field. We clearly observe 
oscillations in the switching time. Moreover there ex- 
ists a region of the parameters (see Fig. [TU|) where the 
switching required more than 16 ns. To understand the 
phenomenon we first plot in Fig. [Tl]the average My mag- 
netisation component as a function of time. The small- 
period magnetisation oscillations are related to the ac- 
ficld oscillations. During the switching process we can 
distinguish two main processes: The first one (occurring 
in the first 2 ns) is related to the efficiency of the magneti- 
sation nucleation in the system while the second one is re- 
lated to the magnetisation relaxation. Fig. [TT] shows that 
at several frequencies the magnetisation reversal process 
started faster in its nucleation part but proceeded slower 
in the relaxation part. 

The nucleation process is associated to efficient spin 
wave generation and then to spinwave instability phe- 
nomena which is a precursor of the magnetisation 
reversal22i^ii£. The role of the ac-field is to excite a 
spinwave with the external frequency. Because of this, 
the data of Fig. [5] resembles the spinwave dispersion rela- 
tion and the ripple configurations - the spinwave modes 
reported in Ref. However, the excited mode may 
be unstable for a given condition and, thus, produce in- 
stabilities and magnetisation reversal. For example, we 
have seen that at these values of the applied field the 
main FMR mode is completely unstable. When we tried 
to excite the homogeneous mode, with the magnetisation 
anti-parallel to the applied field direction, the energy was 
immediately redistributed into inhomogeneous spinwaves 
provoking later the magnetisation reversal. The excita- 
tion of these spinwaves was acting as an additional damp- 
ing source leading to the magnetisation reversal similar 
to the case of Ref. HH. In this case, the excitation with 
the frequencies close to FMR (around 3GHz in Fig. [9]) 
is the most efficient one. The first stage (several ns) is 
also characterized by initial expansion of the domain size, 
especially in the center, as seen in Fig. \7\ 
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FIG. 9: Switching time of the magnetisation as a function of 
the mw frequency for various values of the mw-field amplitude 
Hac and Hac = -31.410e 




v(GHz) 



FIG. 10: Contour plot for the switching time of the ellipsoid 
for applied field Hdc = — 31.410e. 



To illustrate the role of precession, we present in 
Fig. [12] the temporal evolution of the average value of 
the < > component (this value is proportional to the 
energy put into precession by the microwave). The first 
part of these curves characterizes the initial nucleation- 
expansion process (first 2 ns). It is clear that the energy 
is put efficiently into the precession with the frequency 
close to the main FMR mode. This is confirmed by the 
fast magnetisation decay of the < My > component in 
Fig. [TT] and the fast growth of the < 7\f| > value in 
Fig. [T2] This fact is consistent with the hypothesis intro- 
duced in many papers stating that for efficient switch- 
ing the mw frequency should coincide with the FMR 
one^ii^. In the case of magnetic elements behaving as 
one macrospin magnetic moment this would determine 
the overall switching. However, a fast nucleation is insuf- 
ficient for fast magnetisation switching of larger elements. 
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FIG. 11: Temporal evolution for the average < My > mag- 
netisation component (normalized to the saturation value) 
during the microwave-assisted switching process at Hdc = 
-31.41 Oe and Hac = 25.13 Oe. 
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FIG. 12: Temporal evolution for the average < > mag- 
netisation component (normalized to the saturation value) 
during the microwave-assisted switching process at Hdc = 
-31.41 Oe and Hac = 25.13 Oe. 



Fig. [T^l demonstrates that for frequencies close to 6GHz 
(double FMR frequency in the opposite well) the preces- 
sional energy remains constant. We note that only the 
magnetic moments in the centers of domain walls, sepa- 
rating domains with opposite signs (and in the vortex 
centers) are precessing (see Fig. [5]). Thus the expansion 
of domain walls necessary for magnetisation relaxation is 
dependent on the relaxation of these magnetic moments. 
For frequencies and timescales at which the < Af^ > re- 
laxation is very slow, the energy is transferred efhciently 
into the precessional motion and not to the relaxation 
process. Consequently, the further domain expansion is 
extremely slow. We would like to note that this idea 
is in the spirit of the hypothesis introduced in Ref. IsH 
(self-trapping of magnetic oscillations). At the long-time 
scale the magnetisation reversal at this stage proceeds 



again via additional spin-wave generation and reflection 
destabilizing the domain structure. Note also that the 
nucleation process occurs when the magnetisation is anti- 
parallel to the field direction while in the relaxation part 
it is parallel. Thus, the relevant frequencies are different 
in both cases. 

Fig. [12] shows the magnetisation switching time as a 
function of the mw frequency for two dc applied field 
values. As the applied dc field increases, the angle of 
the magnetisation precession increases in the nucleation 
part of the process but decreases in the relaxation part. 
Larger precessional angles lead to nonlinear phenomena. 
Associated with larger precessional angles there is a non- 
linear shift of the frequencies to smaller values^. This 
is in agreement with the results presented in Fig. I13[ 
provided that the relevant frequencies are determined by 
the relaxation process. However, in Fig.lHlthe frequencies 
are shifted to larger values with larger ac-field amplitude. 
This shows that the overall process is much more com- 
plicated and cannot be analysed in terms of the FMR 
frequencies relevant to the switching of one magnetic mo- 
ment only. In fact, different ac-frequencies excite modes 
with different spin wave vector and the instabilities of 
them occur at different threshold amplitudes. 
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FIG. 13: Switching time of the magnetisation as a function 
of the mw frequency for two values of the stationary field Hdc 
and the mw-field amplitude Hac = 25.13 Oe. 



IV. CONCLUSIONS 

Using the micromagnctic model, we have investi- 
gated the mechanism of fast magnetisation switch- 
ing assisted by a linearly polarised microwave field in 
micron-sized magnetic elements. Magnetisation dynam- 
ics in these magnetic patterns is governed by nucleation- 
propagation-relaxation processes and is different to those 
of nanoscale elements behaving as one macrospin. Our 
simulations confirm that the microwave-assisted magneti- 
sation reversal requires a smaller field than in the process 
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under a static applied field. The first contribution to this 
is the deviation of the field from the easy-axis during the 
mw-assisted switching process. However, even when we 
compared the situations when the maximum applied field 
angle is the same in both cases, mw-assisted switching ap- 
peared to require less field. Moreover, our results show 
that at small intensity of the mw-ficld its action is not 
described by the effective field tilt. This happens due 
to the fact that the static and the mw-assisted switching 
processes involve different reversal modes. In the first 
case this mode consists of the nucleation of two vortices 
in opposite corners, their subsequent merging and cre- 
ation of two domain walls which later span the whole 
ellipsoid. In the second case, the reversal mode con- 
sists of a ripple structure, the external frequency being 
responsible for the excitation of a spinwavc mode with 
the corresponding ripple size. For fields below the static 
coercivity field, the magnetisation configuration corre- 
sponds to almost homogeneous background plus excited 
spin wave mode and is metastable. The amplitude of the 
excited spin wave mode grows with applied field and be- 
comes unstable, leading to magnetisation reversal. Thus, 
the spin-wave instability mcchanis m'^^i^^i^° plays a ma- 
jor role for the fast magnetisation switching process in 
nano-size magnetic elements. 

We have shown that the most efficient nucleation is 
at the FMR frequency (and its multiples). However, the 
magnetisation reversal process after the nucleation re- 
quires also an efficient magnetisation relaxation. For ex- 
ample, the domain growth is stimulated by the ac field. 
This happens due to the relaxation of precessing mag- 
netic moments in the domain wall center. In the spirit of 
the mechanism suggested in Ref . [3^. for some frequencies, 
the relaxation process is not efficient, this happens when 
the mw field is coupled to the precessional motion. In this 
case the microwave field efficiently puts the energy into 
precession and not into the propagation-relaxation pro- 
cesses which arc slowed down. As a consequence of the 
interplay of several mechanisms with different relevant 
frequencies, the switching time of magnetic elements is 
a complicated function of the external frequency. The 
above results show that the magnetisation dynamics in 



micron-sized magnetic elements could not be analysed as 
a simple FMR-related phenomenon and neither follows a 
direct correspondence with a spinwave spectrum. 

Finally, we would like to discuss our results in the 
context of available experimental data. Out results are 
in qualitative agreement with experimental observations. 
The Kerr images in Ref. show the creation of ripple 
structures during the mw-assisted switching process in 
a permalloy ellipsoid. Currently there is no data avail- 
able for the switching time of the magnetisation in el- 
lipsoids which could confirm the existence of the oscil- 
lations presented in Fig. [S] However, in several exper- 
imental papers on different magnetic elements we have 
found qualitative similarities with our results. For ex- 
ample, in Co bars the relaxation time observed by the 
fast Kerr technique showed a strong nonmonotonic be- 
havior as a function of applied ficlct. Also in Co bars 
measured by the anisotropic magnetoresistance effecli^, 
the authors have observed the occurrence of several res- 
onance peaks sweeping in frequency. The microwave 
power input necessary for the maximum coercivity re- 
duction was reported to be frequency-dependent in the 
measurements of mw-assisted switching in NiFe mag- 
netic tunneling junctions^. The critical field necessary 
for the microwave assistance was shown to be a compli- 
cated function of the mw frequency with several minima 
for the mw- assisted process going from vortex to onion 
state in magnetic nanorings^^. 
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